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APPENDIX 1 

%%% 

function B = Initcluster(X,m) 

% Get an initial base codebook B at random from input data 
% INPUTS 

% X = input data: each column is a RGB 'vector' 
% m = number of base codevectors 

% OUTPUTS 

% B = base codevector matix 

%need to duplicate some columns of B if X is small 
[n,N] = size(X); 
if(N > m) 

replace = 0; 
else 

replace = 1 ; 
end 

%track what inputs put in B, so no duplication for X large 
chosen = zeros(l ,N); 
B = zeros(n,m); 
for i=l:m 
draw = floor(N*rand +1); 
if(~replace) 
while(chosen(draw)) 

draw = floor(N*rand +1); 
end 
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# 



end 

B(:,i) = X(:,draw); 
chosen(draw) = 1 ; 
end 




function [B,Bq,d,iters,rd,nX] = lgbo(X ? m,base) 

10 % INPUTS 

% X = input data: each column is a RGB 'vector' 

% m = number of codevectors (columns in C) 

% base = number of base codevectors (columns of B) 

15 % OUTPUTS 

% B = base codevector matrix 

% Bq = quantized base codevector matrix 

o/ 0 j = distortion of X when replaced with chosen codevectors 

% iters = # of iterations to reach convergence 

20 

[d,N] = size(X); % d = dimension, N = blocksize 
%random initialization 

Init = initcluster(X 5 base); % choose random initial set 

25 

if base=2 % 1 D (linear interpolation) 

W = [l:(-l/(m-l)):0; 0:(l/(m-l)):l]; %weight matrix 
elseif base=3 %2D 

W=[l 5 0 ? 0;0 ? l 5 0;.25 v 25,.5;-.25,-.25,L5; ... 
30 1,-U;-1,U;0,-1,2;-1A2;]'; 
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elseif base===4 %3D 

W=[eye(4),.25*[2,2J r l;2,-l,24;l,2 r l > 2;-l,l,2,2]]; 
%W=[eye(4) > J25*[3 J 2 f 2J;2,3 f 2J;3 > l,2 f 2;l > 2,3,2]]; 
end 

5 

%variable initialization 
stoppingeps = l.e-5; 
vi - ones(l,m); 

index=zeros(l ,N);mind=index; 
10 cumdist = Lif; lastdist = 0; 

C=Init*W; %interpolate codevectors 



# iterate until convergence 
15 while(abs(cumdist - lastdist) > stoppingeps) 
lastdist = cumdist; 
cumdist = 0; 
iters=iters+l; 



20 while(l) %1 iteration that repeats if B goes singular 
% step (A) 

% form Voronoi regions: for each input, 

% determine which centroid it is closest to 
for i=l:N 
25 V=X(:,i(vi))-C; 

nm=sum(V.*V); %Euclidean distance squared (MSE) 
%input i ! s closest codevector 
[mind(i) ,index(i)] =min(nm) ; 
end 

30 cumdist = sum(mind); 
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# 

f 



%find diagonal matrix N 
for j=l:m 

n(j) = sum(index=j); 

5 end 

% check if B is singular 
% and force it to be non-singular 
num=sum(n~=0) ; 
10 ifnum<base % fewer than base non-zero! 
nz = find(n); 

if (base-num)=l %base=2 
[jyji]=max(mind); 
Init(:,l :(base-l))=B(:,nz); 
15 Init(:,base)=X(:ji); 
else 

[jyji]=sort(mind); 
rnit(:>l :num)=B(: 5 nz); 
jl=N; 

20 for jk=(num+ 1 ) :base 

Init(:jk)=X(:ji01)); 

while(sum(abs(X(:ji(jl))-X(:jiGl-l))))=0) 

jl=jl-l; 
end 

25 j l=j 1- 1 ; 

end % for jk 
end % else 

C=Init*W; %start all over if hit a singular matrix 

else 

30 break; 
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end % if not singular then end while loop 
end%while(l) 

s = zeros(d,m); %get sum matrix S 
5 forj=l:m 

ifn(j)=l 

s(:j) = X(:,index==j); 
elseif n(j) 

s(:j)=sum(X(: 3 index==j) , ) , ; 
10 end 

end %for j=l:m 

Init = s*W'*inv(W*diag(n)* W); %new base codevector matrix B 
C=Init* W; %new codevector matrix C 

15 

end %% while not converged 

nX = B(: 5 index); % save re-constructed block 
d = cumdist; % total block distortion 

20 end % function 
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